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ABSTRACT 

The recent determination of the baryon-to-photon ratio from WMAP data by Spergel 
et al. (2003) allows one to fix with miprcccdentcd precision the primordial abundances 
of the light elements D, ^He, ''He and ^Li in the framework of the standard model 
of big bang nucleosynthesis. We adopt these primordial abundances and discuss the 
implications for Galactic chemical evolution, stellar evolution and nucleosynthesis of 
the light elements. The model predictions on D, '^He and "^He are in excellent agreement 
with the available data, while a significant depletion of ^Li in low-metallicity stars is 
required to reproduce the Spite plateau. 

Key words: Galaxy: abundances - Galaxy: evolution - nuclear reactions, nucleosyn- 
thesis, abundances. 



1 INTRODUCTION 

The standard theory of big bang nucleosynthesis (SBBN; 
e.g., Boesgaard & Steigman 1985; Steigman 1989; Walker 
et al. 1991) accurately predicts the primordial abundances 
of the light elements D, ^He, *He and ''Li, as a function of 
the cosmic baryon density, pb cx: Slbft^, or, equivalently, of 
the baryon-to-photon ratio, r]. Since the baryon density is 
the sole parameter in the SBBN, observations of D, ^He, 
*He and ^Li in astrophysical environments not yet affected 
by subsequent stellar evolution offer a direct way to infer 
the baryon density of the universe and to assess whether 
the SBBN theory correctly describes the first three minutes 
of the hot early universe. Nowadays, there is still disagree- 
ment on which is the value of rj as inferred from i) limits 
on the primordial deuterium abundance from high-redshift 
absorption systems (mostly Damped Lyman a or Lyman 
limit absorbers; e.g., Carswell et al. 1994; Songaila et al. 
1994; Tytler, Fan & Buries 1996; Rugers & Hogan 1996; 
Songaila, Wampler & Cowie 1997; Webb et al. 1997; Buries 
& Tytler 1998; Levshakov, Kegel & Takahara 1998; O'Meara 
et al. 2001; Pettini & Bowen 2001; Levshakov et al. 2002); 
a) estimates of the primordial *He abundance from extra- 
galactic HII regions (e.g., Olive, Skillman & Steigman 1997; 
Izotov et al. 1999; Peimbert, Peimbert & Luridiana 2002; 
see also Pagel 2000 for a review) and in) estimates of the 
pristine ^Li content in old, metal-poor ([Fe/H] < —1.5 dex), 
warm (Teff > 5700 K) dwarf stars in the solar neighbourhood 
(e.g., Sphe & Sphe 1982; Spite, Maillard & Spite 1984; Sphe 



& Sphe 1986; Rebolo, Molaro & Beckman 1988; Thorburn 
1992, 1994; Bonifacio & Molaro 1997; Vauclair & Charbon- 
nel 1995, 1998; Theado & Vauclair 2001; Pinsonneauh et al. 
1999, 2002). Measurements of ^He/H in HII regions in the 
outer parts of the Galactic disc should provide values very 
close to the primordial one, owing to the slow evolution of 
the disc in these regions, but they appear to poorly con- 
strain the 77 range (e.g.. Rood et al. 1998; Bania, Rood & 
Balser 2002). 

Data on cosmic microwave background (CMB) 
anisotropics provide an alternative, independent method for 
constraining 77. The first release of results from the Wilkinson 
Microwave Anisotropy Probe {WMAP; Bennett et al. 2003; 
Spergel et al. 2003) makes the CMB the prime cosmic bary- 
ometer, owing to the high WMAP precision. A combination 
of WMAP data with other finer scale CMB experiments 
(ACBAR - Kuo et al. 2002 - and CBI - Pearson et al. 2002) 
and with astronomical measurements of the power spectrum 
(2dF Galaxy Redshift Survey measurements - Percival et 
al. 2001 - and Lyman a forest data - Croft et al. 2002; 
Gnedin & Hamihon 2002) gives Sib = 0.0224 ± 0.0009, 
or, equivalently, 7710, cmb = 6.1^q' I (Spergel et al. 2003) t, 
where rjio = 10^" r;. By adopting the SBBN predictions and 
using this 77 value, one obtains the primordial abundances 
of D, '^He, *He and ^Li, and can use them to compare with 
observations in low-metallicity environments and to gain in- 
sight into chemical evolution. 

The aim of this paper is to analyse the effect of adopt- 
ing primordial abundances predicted by SBBN constrained 
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Figure 1. Solid lines: primordial abundances of '*He (mass frac- 
tion), D, ^He, ^Li (number ratios relative to hydrogen) as a func- 
tion of )7io = 10^" rj from SBBN. The dashed lines are Icr devi- 
ations. Theoretical predictions are from Fiorentini et al. (1998). 
The data quoted in the Introduction are shown with Icr boxes. 
The vertical band represents the WMAP rjio range. 



by WMAP on the chemical evolution of the Milky Way and 
to gain insights into the mechanisms of hght element pro- 
duction/destruction inside stars. 



2 PRIMORDIAL ABUNDANCES 

In Fig. 1 we show the primordial abundances of D, ^He, 
''He and ^Li predicted by the SBBN model as a function 
of ryio, with their la deviations. Theoretical predictions are 
from Fiorentini et al. (1998). The values inferred from the 
observations quoted in the Introduction are also displayed 
with la boxes. 

The primordial abundance by mass of ''He (Fp) is de- 
rived by analysing the helium content in low-metallicity ex- 
tragalactic Hll regions. Since even in the lowest metallicity 
HII regions some *He must be produced by stars along with 
the heavy elements, it is necessary to correlate Vobs with 
metallicity and to extrapolate to zero metallicity in order 
to infer Yp (Peimbert & Torres- Peimbert 1974). Olive et al. 



(1997) found Yp = 0.234 ± 0.002 for their full data set and 
an even lower value, Yp — 0.230 ± 0.003, for a subset of 
Hll regions with the lowest metallicity. Izotov et al. (1999) 
used high-quality spectroscopic observations of the two most 
metal-deficient blue compact galaxies known, IZwl8 and 
SBS 0335-052, to determine Yp. The weighted mean helium 
mass fraction in the two galaxies, Y — 0.2462 ±0.0015, after 
correction for the stellar ''He production implies a primordial 
''He mass fraction Yp = 0.2452 ± 0.0015. This value corre- 
sponds to rjio — 4.7^0 3; while Olive et al.'s value translates 
into rjio = 1.8 ±0.3. This discrepancy could be less severe, if 
corrections for unseen neutral hydrogen and/or helium ap- 
ply to Izotov et al.'s estimate of Yp (Gruenwald, Steigman & 
Viegas 2002). Very recently, Peimbert et al. (2002) have ob- 
tained Yp = 0.2384 ± 0.0025 from the helium abundances of 
five objects including the three low-metallicity objects with 
the best line determinations and the two objects with the 
lowest metallicity. They used the temperature derived from 
the He 1 lines, which is always smaller than that derived from 
the [O III] lines, and considered the coUisional contribution 
to the Balmer-line intensities, never taken into account in 
the other recent Yp determinations in the literature. Indeed, 
the temperature structure, the coUisional excitation of the 
hydrogen lines, and the ionization structure are problems 
affecting the Yp determination that need to be further anal- 
ysed. The determination of Yp by Peimbert et al. favors rj 
values in the range (2 — 3) x 10"^". 

An alternative approach to derive Yp from Galactic 
globular clusters was suggested by Iben (1968) and Iben & 
Faulkner (1968). This method originally provided low values 
of Yp = 0.23 ± 0.02 (Buzzoni et al. 1983), but was affected 
by significant uncertainties both on the star counts and on 
the stellar evolution theory. Very recent application of the 
method to much more complete data sets and with updated 
stellar models (Cassisi, Salaris & Irwin 2003) indicates in- 
stead Yp ~ 0.243 - 0.244. 

Observations of D in high-redshift QSO absorbers pro- 
vide an estimate of the primordial deuterium abundance. 
The controversy on whether the primordial abundance of 
deuterium was a few times 10~^ or one order of magnitude 
higher (e.g., Carswell et al. 1994; Songaila et al. 1994; Tytler 
et al. 1996; Rugers & Hogan 1996; Songaila et al. 1997; Webb 
et al. 1997; Buries & Tytler 1998) seems now largely solved: 
several low D/H values have been reported in the last few 
years (e.g., Levshakov et al. 1998; O'Meara et al. 2001; Pet- 
tini & Bowen 2001; Levshakov et al. 2002 - Table 1) and 
it has been argued that D/H values higher than a few 10"^ 
are likely to be due to contamination by HI interlopers at 
velocities similar to the isotope shift (e.g., Kirkman et al. 
2001). 

■^He is not a very good baryometer. Rood et al. (1998) 
suggested ('^He/H)p — 1.5lQ g x 10~^ as a reasonable esti- 
mate for the primordial '^He. A smaller range, ('^He/H)p = 
(1.1±0.2) X 10"^ has been more recently obtained by Bania 
et al. (2002) for a Galactic Hll region located at Rg = 16.9 
kpc, namely, for a region where stellar activity should not 
have contaminated the pristine '^He abundance significantly. 
As can be seen in Fig. 1, even this latter, more precise de- 
termination of (^He/H)p does not significantly constrain the 
value of T^io, because of the weak dependence of (''He/H)p 
on the baryon density. However, it points against the low 
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Table 1. D/H measurements at high redshift. 



Authors 



10^ D/H Notes 



Levshakov et al. (1998) 
O'Meara et al. (2001) 

Pcttini & Bowcn (2001) 
Levshakov et al. (2002) 



4.1 — 4.6 Range common to different analyses 

3.0 ± 0.4 Best value from measurements in different objects 

2.54 ± 0.23 Lyman limit system at z^bs = 2.536 along the line of sight of HS 0105+1619 

1.65 ± 0.35 DLA at Zabs = 2.0762 along the Une of sight of Q 2206-199 

3.75 ± 0.25 DLA at z^bs = 3.025 along the line of sight of Q 0347-3819 



7710 range (?7io < 3) suggested by the high-D, low-'' He deter- 
minations. 

The actual value of the primordial ^Li abundance is still 
controversial, with some people favoring the hypothesis that 
the so-called Spite plateau is truly representative of the pris- 
tine ^Li content of halo dwarfs, and hence of the primordial 
■^Li abundance (Spite & Spite 1982; Spite et al. 1984, Spite & 
Spite 1986; Rebolo et al. 1988; Spite et al. 1996; Bonifacio & 
Molaro 1997), and other authors suggesting the existence of 
non-standard depletion processes to reproduce the observed 
plateau value starting from a higher primordial one (e.g., 
Thorburn & Beers 1993; Vauclair & Charbonnel 1995, 1998; 
Theado & Vauclair 2001; Pinsonneault et al. 1999, 2002). 
On the other hand, the primordial lithium abundance could 
be even lower than that shared by warm halo dwarfs because 
of significant early contribution from Galactic chemical evo- 
lution (a + a fusion mechanism and stellar nucleosynthesis; 
e.g., Ryan, Norris & Beers 1999; Suzuki, Yoshii & Beers 
2000). A set of representative primordial abundances of ^Li 
is shown in Fig. 1. Notice that the most recent estimate 
of the plateau value from Bonifacio (2002), loge(^Li)p'i" = 
2.317 ± 0.014io- ± O.OSsys, is not shown, because it refers 
to a preliminary data analysis published inso far only as 
conference proceedings. The quoted value results from the 
analysis of a high-quality sample of 22 stars satisfying the 
conditions Tes > 5700 K, [Fe/H] < —1.5, taking into ac- 
count the effects of standard depletion and non-LTE. This 
value is higher than previous estimates by Spite & Spite 
(1982) [loge(^Li)p = 2.05 ± 0.15] and Bonifacio & Molaro 
(1997) [loge(^Li)p = 2.238 ± 0.012i<, ± 0.05sys], but always 
lower than those suggested if non-standard stellar depletion 
is acting in low-mass, warni halo stars. 

The first release of results from WMAP allows us to 
independently constrain the range of values of r)io with un- 
precedented precision. A value of 7710 = 6.11° 2 is found 
(Spergel et al. 2003). The rjio range is shown in Fig. 1 as a 
vertical band. It is clearly seen that the low-^He and high- 
D values are completely ruled out. In particular, for D this 
was already a prediction from Galactic chemical evolution 
(GCE) models (Steigman & Tosi 1992; Galh et al. 1995; 
Prantzos 1996; Tosi et al. 1998; Chiappini, Renda & Mat- 
teucci 2002). In fact, the level of stellar astration required 
by chemical evolution models for deuterium to match its 
present interstellar medium (ISM) abundance starting from 
a high primordial value is inconsistent with a number of 
observational constraints for the Milky Way. 



t log£(^Li) = log(N7Li/NH) + 12, where Ntl; and Nh axe the 
abundances by number of ^Li and H, respectively. 



In the following, we adopt the primordial abundances 
resulting from the 7710 range derived from WMAP data and 
study the evolution of the light elements D, ^He, ''He and ^Li 
in the Galaxy. We discuss the implications of these model 
results on stellar evolution and nucleosynthesis. 



3 LIGHT ELEMENT EVOLUTION IN THE 
MILKY WAY 

We adopt two different chemical evolution models for the 
Milky Way, the one developed by Chiappini et al. (1997, 
2001) and that developed by Tosi (1988a, b - model called 
Tosi 1). Both models satisfy the main observational con- 
straints for the solar neighbourhood as well as the whole 
Galactic disc. In the following, we underline the major dif- 
ferences and similarities between the models. 

In both models, the Galactic disc is divided into con- 
centric rings without radial flows between them and the disc 
of the Galaxy is built up inside-out (Mattcucci & Francois 
1989) from gas infalling from outside (see also Larson 1976; 
Tinsley & Larson 1978; Tosi 1982). In Chiappini et al.'s 
model, the infall rate is exponentially decreasing in time, 
with an e-folding time increasing with increasing Galactic 
radius. In particular, td(-Rg) = 1.03 x Rg — 1.27 Gyr (Ro- 
mano et al. 2000) . This leads to a faster disc formation in the 
innermost regions. The infalling gas has a primordial chem- 
ical composition. In the model of Tosi, a constant (in time) 
and uniform (in space) infall density of 0.004 Mq kpc~^ 
yr^'^ is assumed, so that more mass is accreted in the outer 
disc than in the inner one; the infall law is formally expo- 
nential and decreasing in time, but practically constant. The 
infalling gas is slightly metal-enriched, Zinf = 0.2 Zq (see 
Tosi 1988b). In Chiappini et al.'s model, the star formation 
rate at a given radius Rq and time t is explicitly depen- 
dent on Egas(iiG,i), Etot(J?G,i) and Etot(J?o, toai), i.e., the 
gas and total mass densities at that radius and time and the 
present-day total mass density at the solar position. In Tosi's 
model the star formation is explicitly dependent on the gas 
and total mass density currently observed at each Galacto- 
centric distance and exponentially decreasing in time (with 
e-folding time 15 Gyr). In the model of Chiappini et al., a 
threshold in the star formation process is considered, which 
makes the star formation in the disc go to zero when the gas 
density falls below a critical value (S**^ = 7 Mq pc~^). This 
threshold is not considered in Tosi's model. The two mod- 
els differ also in the adopted initial mass functions (IMFs): 
Chiappini et al. adopt Scale's (1986) IMF, while Tosi adopts 
Tinsley's (1980) IMF. Both models avoid the instantaneous 
recycling approximation, i.e., they account for the stellar 
lifetimes. However, they adopt different stellar lifetimes. 
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Figure 2. Evolution of D/H in the solar neighbourhood. The 
thick continuous line are the predictions of Model CRM; the thin 
continuous line are the predictions of Model Tosi 1 (see text). 
The vertical bars at t = represent the range of primordial D 
abundance from WMAP (smaller bar) and high-redshift QSO ab- 
sorbers (bigger bar - Levshakov et al. 1998; O'Meara et al. 2001; 
Pettini & Bowen 2001; Levshakov et al. 2002). The PSC data 
(Geiss & Gloeckler 1998) and the local ISM value (Linsky 1998) 
are also shown, at 1 and 2a (thick and thin bars at t = 8.5 and 
13 Gyr, respectively). The vertical dotted line at t = 13 Gyr rep- 
resents the range of abundances derived from Copernicus, HST- 
GHRS, IMAPS, STIS and FUSE data (see, e.g., Vidal-Madjar 
et al. 1998; Jenkins et al. 1999; Moos et al. 2002; Hoopes et al. 
2003). 



For the stellar nucleosynthesis of D, "^He and *He we 
assume that D is only destroyed in stellar interiors^, while 
^He and *He are partly produced and partly burnt to form 
heavier species (e.g., Dearborn, Steigman & Tosi 1996). The 
evolution of '^He is assumed to be strongly affected by extra- 
mixing and subsequent Cool Bottom Processing (CBP) in 
low-mass stars as prescribed by Sackmann & Boothroyd 
(1999a, b). These prescriptions are the same as those already 
adopted by Galh et al. (1997), Palla et al. (2000) and Chiap- 
pini et al. (2002) to reproduce the observed abundances of D 
and ■^He. For ^Li, the adopted nucleosynthesis prescriptions 
are those described by Romano et al. (1999, 2001). 



3.1 



Evolution of D, ^He and ''He in the Galaxy 



Our models assume (D/H)p = 2.5 x 10"^ (^He/H)p = 
0.9 X 10~* and = 0.248, consistently with the range of -q 
inferred from WMAP data. These adopted values happen to 
be the same as adopted by Chiappini et al. (2002) for their 
Model C-II, except for ^He, for which we assume a slightly 
lower primordial abundance. Model C-II also assumes that 
93 per cent of the stars with mass M < 2 Mq suffer some 
extra-mixing, resulting in an overall "^He destruction (see 
Chiappini et al. 2002 for details; see also Galli et al. 1997; 
Charbonnel & do Nascimento 1998). We assume the same 



5i A mechanism for post big bang deuterium production (interac- 
tion of protons accelerated in flares with the stellar atmosphere 
creating a flux of free neutrons, which undergo radiative capture 
on atmospheric protons resulting in D synthesis) has been pro- 
posed (MuUan & Linsky 1999), but we do not include it here. 
Indeed, it has been recently proved that this production channel 
does not allow for D production at a level which will reverse the 
monotonic decline of D (Prodanovic & Fields 2003). 
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Figure 3. Upper panel: evolution of "^He/H in the solar neigh- 
bourhood. Thick continuous line: Model CRM; thin continuous 
line: Model Tosi 1; dotted line: the same as Model Tosi 1, but 
assummg (D/H)p = 20 X 10-^5; dashed line: the same as Model 
Tosi 1, but assuming (D/H)p = 10 X 10~^; dot-dashed line: the 
same as Model Tosi 1, but with 100 per cent of low-mass stars 
suffering extra-mixing in conjunction with CBP occurrence. Data 
(vertical bars at t = 8.5 and 13 Gyr) are from Geiss & Gloeckler 
(1998). Lower panel: radial distribution of '^He/H for the same 
models in the upper panel. The dots are HII region data from 
Bania et al. (2002). 



percentage throughout this paper as well, unless otherwise 
stated. 

The adopted theoretical SBBN predictions are from 
Fiorentini et al. (1998). The assumed (^He/H)p value is 
the lowest still compatible with WMAP data (at la level 
- see Fig. 1). This choice is motivated by the fact that 
GCE models always tend to overproduce ^He (especially 
if extra-mixing and CBP are not included in the compu- 
tations - see, for instance, figure 4 of Palla et al. 2000). 
Therefore, it is meaningful to start from the lowest ad- 
missible value for (^He/H)p. It is worth noticing that dif- 
ferent SBBN computations do not always lead to val- 
ues of (D/H)p, (^He/H)p, Fp and (''Li/H)p which agree 
(within the errors) with those adopted here. For instance, by 
adopting the SBBN predictions by Cyburt, Fields & Olive 
(2001), in light of the WMAP determination of r) one gets 



(D/H)p = 2.7AtlS, X 10-^ (^He/H)p 



Fp = 0.2484+»-°°o^' 



9.30tJ;^? X 10" 



p - ...,o,_o.ooo5 and ('Li/H)p = ^i.lGt'olt x 10"^". In- 
stead, for the OSU code, at 7)w = 6.14, (D/H)p = 2.57x 10"^ 
and (^Li/H)p = 4.51 x 10~^° (G. Steigman, private commu- 
nication) . These results are fairly different from those listed 
above, but similar to the ones from the Buries et al. (1999) 
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code [(D/H)p = 2.60 x 10"^ and (■^Li/H)p = 4.91 x lO"^"]. 
However, we have checked that no appreciable changes in 
the GCE predictions are produced by adopting these values 
rather than those used here. 

In Fig. 2 we show the temporal behaviour of the deu- 
terium abundance in the solar neighbourhood predicted by 
Model CRM (the same as Model C-II of Chiappini et al. 
2002, but with a lower ^He primordial abundance; thick 
continuous line) compared to that obtained by Model Tosi 
1 (thin continuous line) starting from the same primordial 
deuterium abundance and considering extra-mixing in 93 
per cent of low-mass stars. It can be seen that, notwith- 
standing the differences between the two models mentioned 
above, the predicted behaviour is almost the same. In par- 
ticular, the two models predict nearly the same depletion 
factor for deuterium, ^D,p/^D,tGai ^ 1-5, where Xo.p and 
Xd, toai 8.re the primordial and present-day D abundances 
by mass, respectively. The vertical bars at f = represent 
the range allowed for the primordial deuterium abundance 
from WMAP data (smaller bar) and from high-redshift QSO 
absorbers (bigger bar - Levshakov et al. 1998; O'Meara et 
al. 2001; Pettini & Bowen 2001; Levshakov et al. 2002). The 
vertical solid bars at t — 8.5 Gyr are the Geiss & Gloeck- 
ler (1998) Protosolar Cloud (PSC) data at 1 and 2a (thick 
and thin lines, respectively), while the vertical solid bars at 
t = 13 Gyr represent the local ISM value by Linsky (1998; at 
1 and 2a - thick and thin solid lines, respectively) . The dot- 
ted line at t = 13 Gyr shows the most likely range of D/H 
variation as allowed from measurements with Copernicus, 
HST-GHRS, IMAPS, STIS and FUSE along several lines 
of sight (e.g., Vidal-Madjar et al. 1998; Jenkins et al. 1999; 
Moos et al. 2002; Hoopes et al. 2003). However, notice that 
values as high as D/H = 4 x 10""" have also been reported 
from HIE observations (see Vidal-Madjar et al. 1998 for ref- 
erences, but see also Vidal-Madjar & Ferlet 2002). Starting 
from the WMAP primordial D determination, the protoso- 
lar and local deuterium abundances are very well reproduced 
by assuming that D is only destroyed in stars. 

In Fig. 3, upper panel, the evolution of ''He/H in the so- 
lar neighbourhood predicted by the same models of Fig. 2 is 
displayed, together with predictions from other three mod- 
els, which differ from Model Tosi 1 either in the assumed pri- 
mordial deuterium abundance or in the percentage of low- 
mass stars suffering extra-mixing. In particular, the cases 
(D/H)p = 20 X 10-^ (D/H)p = 10 X 10"^ and 100 per cent 
of low-mass stars destroying '^He through CBP are shown, 
as a dotted, dashed and dot-dashed thin line, respectively. 
The protosolar and local ^He abundances are from Geiss & 
Gloeckler (1998). It is seen that models starting from a high 
primordial deuterium abundance are not able to reproduce 
the protosolar and/or local ^He abundance, even allowing 
for ~ 90 per cent of extra-mixing in low-mass stars. This 
gives further support to the argument that the primordial D 
abundance cannot be higher than (D/H)p ~ 4 x 10~^ (Chi- 
appini et al. 2002; see also Tosi et al. 1998). A model where 
100 per cent of low-mass stars destroy their ^He through 
CBP results in an almost flat ^He behaviour from the time 
of the protosolar nebula formation to now. This is clearly an 
extreme case, since the high ^He abundances derived for a 
small sample of Galactic planetary nebulae (e.g., Balser et 
al. 1997; Balser, Rood & Bania 1999; Palla et al. 2000) sug- 
gest that at least a minor fraction of the stars of 1 - 2 Mq do 
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Figure 4. Temporal (upper panel) and spatial (lower panel) vari- 
ation of the helium isotopic ratio in the solar neighbourhood and 
across the Galactic disc at the present time, respectively. The 
predictions from Models CRM (thick solid lines) and Tosi 1 with 
93 and 100 per cent of extra-mixing (thin solid and dot-dashed 
lines, respectively) are compared to the observations (see text for 
details). 



indeed produce ^He in significant quantities. In Fig. 3, lower 
panel, the behaviour of ^He/H across the Galactic disc as 
predicted by the same models for the present epoch is dis- 
played and compared to HII region observations by Bania 
et al. (2002). We do not recover the flat behaviour suggested 
by Bania et al. (2002), but rather predict a negative ^He gra- 
dient in the 4-18 kpc Galactocentric distance range, unless 
the case of 100 per cent of low-mass stars suffering extra- 
mixing is considered. 

As far as ^He is concerned, starting from the WMAP 
value of Yp = 0.248, a value of Yq = 0.261 (Model CRM) 
and Yq = 0.27 (Model Tosi 1) is found at the time of Sun 
formation. The latter nicely compares to the value at birth 
of the Sun oiY = 0.275±0.01 from Grevesse & Sauval (1998; 
see also Bahcall, Pinsonneault, & Basu 2001). Notice that 
by adopting the Tinsley (1980) IMF, i.e. the same IMF as- 
sumed by Model Tosi 1, Model CRM predicts Yq = 0.274. A 
detailed discussion of the effects of changing the IMF in the 
chemical evolution code will be presented in a forthcoming 
paper. 

In Fig. 4 we show the temporal (upper panel) and spa- 
tial (lower panel) variation of the helium isotopic ratio in 
the solar neighbourhood and across the Galactic disc at the 
present time, respectively. Only the predictions from Mod- 
els CRM and Tosi 1 with 93 per cent and 100 per cent of 
extra-mixing are shown. The data for the solar nebula and 
the local ISM (vertical bars in the upper panel) are taken 
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Figure 5. a) Evolution of ^Li in tlie solar neiglibourliood. Tlie continuous line is a model adopting the prescriptions of the best model 
of Romano et al. (2001), except for the fact that it starts from a higher primordial lithium abundance of log£('^Li)p ~ 2.6. Data are 
non-LTE values for dwarf stars hotter than 5700 K from Bonifacio & Molaro (1997), Chen et al. (2001) and Ryan et al. (2001). The solar 
(Grevesse & Sauval 1998), meteoritic (Nichiporuk & Moore 1974) and TTauri star (at 2cf error; Stout-Batalha et al. 2000) abundances 
are also shown with different symbols, b) Evolution of ^Li in the solar neighbourhood in a log£(^Li) versus age diagram. The chemical 
evolution model is the same as in Fig. 5a. Data for Galactic open clusters and globular clusters (big circles) are shown with their error 
bars as well as different measurements of the ISM value (triangles) (see text for references). The meteoritic value (pentagon) and data 
for field dwarfs with [Fe/H] > — 1.4 (small filled circles) are also shown. 



from Geiss & Gloeckler (1998; at 1 and 2cr - thick and thin 
lines, respectively) and Salerno et al. (2003; Icr-bar on the 
right at t = 13 Gyr). The measurement of '^He/''He in the 
local ISM from the COLLISA experiment on board of MIR 
[='He/''He = (1.7 ± 0.8) x 10"**; Salerno et al. 2003] in- 
volves neutral gas reaching Earth's orbit before interacting 
with the solar EUV photons, which thus keeps the origi- 
nal isotopic abundance ratios. If heavier weight is given to 
this more recent measurement, a case where nearly all low- 
mass stars suffer extra-mixing should be preferred. A reli- 
able estimate of the '^He/'*He gradient across the Galactic 
disc is hence highly desirable in order to establish which is 
the actual fraction of low-mass stars destroying their pristine 
^He. In fact, different models, assuming different degrees of 
extra-mixing, predict large differences in the ^He/^He ra- 
tio across the Galactic disc, especially at inner radii (Fig. 4, 
lower panel). Assessing this point could have important con- 
sequences on our knowledge of the mechanisms of ^He pro- 
duction/destruction inside stars. Unfortunately, the value of 
the ^He/^He ratio at different positions along the Galactic 
disc determined by observations is very uncertain, owing to 
the uncertain, model-dependent ionization corrections which 
apply to the raw data. 

The temporal variation of D/H and "^He/H in the so- 
lar neighbourhood predicted by the two different chemical 
evolution models (Model CRM and Model Tosi 1), when 
starting from the same primordial abundances and adopt- 
ing the same nucleosynthesis prescriptions, is fairly simi- 
lar (see Figs. 2 and 3). However, Model Tosi 1 predicts a 
slightly lower deuterium abundance and a slightly higher 



''He abundance from t ~ 5 Gyr to now, except for the last 
1 Gyr of evolution, when a higher '^He mass fraction in the 
ISM is expected according to Model CRM than according to 
Model Tosi 1. The same behaviour is recovered also in the 
^He/''He versus time plot (Fig. 4). The differences in the 
abundance trends can be explained as due to the different 
stellar lifetimes, IMF slopes, infall and star formation laws 
adopted by the two models. In particular, in Model Tosi 1 
stars in the mass range 0.5-0.9 Mq never die, while they 
do in Model CRM, and these low-mass stars contribute an 
important fraction of ^He at late times. The different choice 
of the IMF produces a smaller - although still non-negligible 
- effect on the predicted ^He evolution. 

3.2 Evolution of ^Li in the solar neighbourhood 

Once the SBBN model predictions from Fiorentini et al. 
(1998) are adopted, the WMAP rjio value of 6.1 results in 
a primordial lithium abundance of (^Li/H)p ~ 3.75 x 10~^" 
(see Fig. 1). In the usual loge('^Li) notation, this translates 
to loge(^Li)p ~ 2.6, higher than any previous estimate of 
the primordial ^Li abundance from observations of warm, 
metal-poor halo dwarfs under the hypothesis that they have 
preserved all their pristine lithium and have not been signif- 
icantly polluted by early Galactic lithium production. In the 
following, we analyse the effect of adopting such a high value 
for the primordial ^Li abundance in the chemical evolution 
code. 

In Figs. 5a and 6a, we show the log£('^Li) versus [Fe/H] 
diagram for field dwarfs in the solar neighbourhood hot- 
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Figure 6. Same as Fig. 5, but here we show also the single contributions to the Li enrichment by the different stellar lithium factories 
we consider: AGB stars (short-dashed line); SNell (dotted line); low- mass red giants (long-dashed line); novae (dot-dashed line). 



ter than 5700 K (data are from Bonifacio & Molaro 1997; 
Ryan et al. 2001; Chen et al. 2001). The solar (Grevesse & 
Sauval 1998), meteoritic (Nichiporuk & Moore 1974) and 
TTauri (Stout-Batalha, Batalha & Basri 2000) values are 
shown as well. In the context of standard stellar evolu- 
tion, warm halo stars ([Fe/H] < —1.5 dex) should not dis- 
play lithium abundances significantly lower than that in 
the gas out of which they formed (Deliyannis, Domarque 
& Kawaler 1990). In fact, in standard stellar evolution mod- 
els lithium is destroyed mainly during the pre-main sequence 
phase (during the main sequence lifetime, the convective re- 
gions are too shallow to allow Li depletion) and pre-main 
sequence Li depletion is not expected at low mctallicitics 
(D'Antona & Mazzitelli 1984). However, the halo stars could 
destroy a significant fraction of their initial lithium content, 
if they were suffering for non-standard destruction mecha- 
nisms (e.g., Vauclair & Charbonnel 1995, 1998; Theado & 
Vauclair 2001; Pinsormcault ct al. 1999, 2002). Since these 
non-standard depletion mechanisms are responsible for a 
broad spread in the Li abundances, as actually seen in the 
data at [Fe/H] > —1.0 dex, the remarkable flatness and thin- 
ness of the Spite plateau represent a major challenge to stel- 
lar evolution models taking non-standard ^Li depiction into 
account. Recent models considering gravitational waves as 
the major depletion cause (Talon & Charbonnel 2003) seem 
indeed capable of reproducing both the observed '^Li spread 
at high metallicities and the tightness of the Spite plateau 
at low ones. These models predict that the observed Spite 
plateau is actually 4-5 times lower than the primordial value 
(Charbonnel 2003, private communication). 

In Figs. 5b and 6b, we show the log e('^Li) versus age di- 
agram, which allows for a bettor appreciation of the lithium 
variations in the last 1 Gyr. Data for Galactic open clus- 
ters as well as old globular clusters are shown. The big filled 
circles represent mean lithium values, while the big open 



circles represent the highest lithium abundance observed in 
the cluster. Different measurements of the ISM value are 
also shown. For NGC 6397, the mean value for 12 turnoff 
stars corrected for non-LTE and standard depletion is that 
by Bonifacio et al. (2002). For 47Tuc, the ^Li abundance is 
the mean value for 2 turnoff stars determined by Pasquini & 
Molaro (1997). For M 92, the highest non-LTE Li abundance 
measured by Boesgaard et al. (1998) is shown. The age of all 
these clusters is taken from Carretta et al. (2000). For M 67, 
data are from Pasquini, Randich & Pallavicini (1997); for 
NGC 752, from Balachandran (1995 - average of the three 
stars with the largest Li abundances on the hot side of the 
Li dip). For NGC 3680, wo assume the highest lithium abun- 
dance from Pasquini, Randich & Pallavicini (2001); for the 
Hyades, the mean value is taken from Balachandran (1995 - 
average of the three stars with the largest Li abundances on 
the cool side of the Li dip). The values for the Pleiades and 
a Per arc from Soderblom et al. (1993) and Randich et al. 
(1998), respectively, while those for IC2391 and IC2602 are 
non-LTE values from Randich et al. (2001). For IC4665 and 
NGC 2264, non-LTE values from Martm & Montes (1997) 
and Soderblom et al. (1999), respectively, are displayed. The 
value for a Ori is the moan value from Zapatcro Osorio et al. 
(2002). The ISM values are from Lemoine et al. (1993; Une 
of sight towards p Oph) and Knauth et al. (2000; lines of 
sight towards o Per and ( Per). Data for field dwarfs (Chen 
ct al. 2001) are also shown for comparison. 

The thick continuous lines in Figs. 5 and 6 represent the 
predictions of a model where Galactic cosmic ray (GCR) and 
stellar lithium production are taken into account in order to 
rise the ^Li content in the ISM from its primordial value of 
log e(^Li)p ~ 2.6 to the meteoritic [log e('^Li) — 3.3] and local 
[log£(^Li) ~ 3.2] ones. In particular, lithium is mostly pro- 
duced by low-mass, long-lived stellar sources (low-mass stars 
on the red giant branch and novae). Only a minor contribu- 
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tion comes from stars on the asymptotic giant branch and 
from Type II supernovae (see Romano et aL 1999, 2001 and 
Fig. 6) contrary to what has been suggested by Travaglio 
et al. (2001), who invoke a large contribution from AGB 
stars (see, however, Ventura, D'Antona & Mazzitclli 2002 
for a critical analysis of their results). It is found that the 
amount of stellar production required in order to match 
the meteoritic data does not change with changing the as- 
sumed primordial '^Li abundance from loge(^Li)p ~ 2.2 to 
loge(^Li)p ~ 2.6, owing to the fact that the evolution of ^Li 
during almost the whole Galgixy evolution is practically de- 
termined solely by the amount of lithium produced through 
stellar and GCR processes (cf. figure 6 of Romano et al. 
2001). The different lines in Fig. 6 represent the predictions 
of models where lithium is contributed by only a single stel- 
lar Li factory: asymptotic giant branch (AGB) stars (short- 
dashed line); Type II supernovae (SNell; dotted line); low- 
mass red giants (long-dashed line); novae (dot-dashed line). 

We conclude that, even in the case in which the primor- 
dial abundance of lithium is as high as loge(''Li)p ~ 2.6 (as 
suggested from recent WMAP data), our previous conclu- 
sions on the ^Li evolution in the solar neighbourhood are left 
unchanged; in particular, the rise off the primordial plateau 
value is still explained as due to the same important, late 
'^Li contribution from long-lived stellar sources (low-mass 
red giants and novae; Romano et al. 1999, 2001). However, 
if the loge(^Li)p value is as high as ~ 2.6, we strongly need 
some non-standard depletion mechanism able to explain the 
absence of scatter and the high level of flatness observed in 
halo stars over a quite large range of metallicity. 



4 FINAL REMARKS AND CONCLUSIONS 

The recent results from WMAP and their direct conse- 
quences on the primordial abundances of the light elements 
have finally allowed us to check that GCE models able to 
reproduce all the major observed properties of the Milky 
Way axe consistent also with SBBN predictions. It has been 
shown several times in the last decade (e.g., Stcigman & Tosi 
1992; GaUi et al. 1995; Prantzos 1996; Tosi 1996; Tosi et al. 
1998; Chiappini et al. 2002) that only a very moderate D 
depletion from its primordial abundance to the present one 
is allowed to let the models reproduce the observed radial 
distributions of chemical abundances, star and gas densities 
and star formation rates, as well as the age-metallicity rela- 
tion and the G-dwarf metallicity distribution. However, un- 
til now there was no definitive observational evidence on the 
primordial values, since one cannot completely exclude that 
even high-redshift, low- metallicity QSO absorbers might be 
already polluted by stellar nucleosynthesis. Now we have 
seen that the D, ^He and ''He produced during the big bang 
are in excellent agreement with Galactic evolution require- 
ments. For "^Li, some problems do arise, but they mainly 
concern our understanding of the mechanisms of lithium di- 
lution/destruction in stars rather than chemical evolution. 
Indeed, by assuming the WMAP primordial abundance of 
log c('Li)p ~ 2.6, we still need the same important, late con- 
tribution of Li from long-lived stellar sources (low-mass red 
giants and novae) required in order to explain the observa- 
tions in meteorites when starting from the lower value of 
loge('^Li)p ~ 2.2 suggested from observations of halo stars 



under the hypothesis that they neither destroyed their pris- 
tine lithium nor suffered any pollution by early GCE (see 
also Romano et al. 2001). In fact, the amount of lithium pro- 
duction from stars and GCRs required by our model in order 
to explain the observations is such that it largely ovorwelms 
the lithium primordial abundance. However, if the adopted 
primordial value of log£r('^Li)p ~ 2.6 is confirmed, a mech- 
anism able to deplete Li in halo stars while preserving the 
flatness of the Spite plateau and producing almost no scatter 
in the star to star abundance at low motallicities is strongly 
needed. Gravitational waves seem interesting in this respect 
(see Talon & Charbonnel 2003), but we still wait for full 
metallicity dependent computations. 

We should recall at this point that our standard cficmi- 
cal evolution models refer to large-scale, long-term phenom- 
ena and cannot account for small-scale, short-term varia- 
tions. If we want to reproduce also the observed spread in 
the abundances of the light elements, as observed, for in- 
stance, by FUSE in the local ISM (e.g.. Moos et al. 2002), 
or for ^He/H across the Galactic disc (Bania et al. 2002), we 
should also take into account the possible inhomogeneities 
in the chemical enrichment of each region and the effects of 
possible orbital diffusion of the stars. 
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